On quantum microcanonical equilibrium 



yo". Dorje C. Brody 1 , Daniel W. Hook 2 , and Lane P. Hughston 3 

f"*) ' 1 Department of Mathematics, Imperial College, London SW7 2BZ, UK 

2 Blackett Laboratory, Imperial College, London SW7 2BZ, UK 
CN| ■ 3 Department of Mathematics, King's College London, The Strand, London WC2R 2LS, UK 

o 

' Abstract. A quantum microcanonical postulate is proposed as a basis for the equilibrium 

(-H . properties of small quantum systems. Expressions for the corresponding density of states are 

' derived, and are used to establish the existence of phase transitions for finite quantum systems. 

A grand microcanonical ensemble is introduced, which can be used to obtain new rigorous 
results in quantum statistical mechanics. 
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1. Introduction 

The purpose of this paper is to examine properties of quantum systems in thermal equilibrium. 
Questions that arise in this context, for example, are: "What is the state of a system in 
equilibrium?" or "What is the temperature of an isolated system in equilibrium?" In the 
case of a classical system immersed in a heat bath, the equilibrium distribution takes the Gibbs 
form exp(—f3H)/Z((3), where (3 = 1/ksT is the inverse temperature of the bath. What about a 
quantum system? Is the equilibrium state given by the Gibbs density matrix exp(— /3ff)/Z(/3)? 
If so, how does one verify that the parameter (3 appearing in the density matrix is the inverse 
^ { - temperature of the bath? 

To investigate questions of this kind it is useful to consider first the classical situation. In the 
classical case, we take the system and the bath as a whole and regard this as a single isolated 
^ . system. The Hamiltonian (symplectic) structure of classical phase space T then allows us to 

define the density of states fl(E) = j r 5(H(x) — E)dV as the weighted volume of the phase space 
occupied by states with energy E. In equilibrium the state of the system maximises entropy 
and thus is given by a uniform distribution over the energy surface; this can be derived if the 
Hamiltonian evolution exhibits ergodicity. The entropy of the equilibrium state is thus given by 
S{E) = ks ln£l(E), and the temperature is defined by the thermodynamic relation TdS = dE. 
These are the necessary ingredients for the consideration of the equilibrium properties of a small 
subsystem. In particular, under a set of reasonable assumptions, it is possible to deduce, by use 
of the law of large numbers, that the equilibrium properties of a small subsystem are described 
by the Gibbs state. A complete derivation of these results is outlined in the seminal work of 
Khinchin [lj. Although the derivation of the equilibrium state is surprisingly complicated, once 
the relevant assumptions are specified, there are no ambiguities in the matter, and familiar 
results associated with the canonical ensemble can be obtained rigorously. 

The situation is markedly different in the case of a quantum system. First, in the usual Hilbert 
space formulation of quantum mechanics it is not clear how one can exploit the Hamiltonian 
structure. This leads to a difficulty in defining the temperature of a closed system. Second, 
since no rigorous derivation of the temperature exists (at least for finite quantum systems), it 



is not possible to verify whether the parameter (3 appearing in the Gibbs density matrix agrees 
with the inverse temperature of the bath. 

In the literature on quantum statistics it is often postulated that the microcanonical density 
matrix of a quantum system with eigenenergy Ei is given by the projection operator onto the 
Hilbert subspace spanned by states with that energy, normalised by the dimensionality 
of that subspace. The entropy is then defined by the expression S = ks lnre^. A rigorous 
derivation of this density matrix is given by Khinchin [2j; however, the assumptions required to 
obtain the result go beyond those required for the classical case. In particular, it is necessary to 
forbid all superpositions of states with different energy. The exclusion of general superpositions, 
however, contradicts the superposition principle of quantum mechanics. This incompatibility 
between quantum mechanics and quantum statistical mechanics is an issue that has troubled 
many authors. For example, Schrodinger remarked in this connection that ". . . this assumption 
is irreconcilable with the very foundations of quantum mechanics" , and that ". . . to adopt this 
view is to think along severely 'classical' lines" [3]. Confronted with this apparent contradiction, 
Schrodinger was nonetheless able to offer an argument to show, in effect, that in thermodynamic 
limit (where the number of particles in the system approaches infinity) the assumption that 
general superpositions are forbidden is justified [3]. 

There is another important shortcoming in the familiar derivation of quantum statistical 
mechanics, namely, that the entropy is a discontinuous function of the energy. As a consequence, 
the temperature of a finite isolated system is undefined. This issue is addressed by Griffiths [3J, 
who demonstrated the existence of a thermodynamic limit in which thermodynamic functions 
are well defined. We thus see that to make sense of the conventional approach to quantum 
statistics, a "macroscopic" limit is required. In this limit, however, we expect quantum systems 
to behave semiclassically so that superpositions, in particular, are excluded. For finite quantum 
systems, these issues remain unresolved. 

While the notion of a thermodynamic limit was justified both theoretically and experimentally 
some forty years ago, there have been experiments carried out on quantum systems over the past 
decade that involve small numbers of particles (see Gross [5] and references cited therein). In 
particular, phase transitions have been observed in small systems — for example, the spherically 
symmetric cluster of 139 sodium atoms exhibits a solid-to-liquid phase transition at about 
267 K [6j. Such experiments demonstrate the breakdown of the conventional approach in which 
phase transitions are predicted only in thermodynamic limits. 

To obtain an equilibrium distribution that is well defined for finite systems, and to address 
the issue of the observed finite-size phase transitions, we have recently introduced an alternative 
formulation to quantum microcanonical equilibrium [TJ. The idea is to follow the derivation of 
the traditional result, as outlined in Khinchin 2j, as closely as possible, but to relax just one of 
the assumptions; namely, for a fixed energy E, we allow the system to be in a superposition of 
energy eigenstates with distinct eigenvalues. 

2. Thermodynamic equilibrium 

The idea of the new microcanonical equilibrium can be described heuristically as follows. We 
consider a gas consisting of a large number N of weakly-interacting identical quantum molecules. 
As in the conventional approach, the intermolecular interactions are assumed strong enough to 
allow the gas to thermalise but weak enough so that, to a good approximation, the total system 
energy can be written as Y2i=i^i ~ ^totab where {Hi\i=i,2,...,N ar e the Hamiltonians of the 
individual constituents. If the composite system is in isolation, then the total energy is a fixed 
constant: Yli=i(Hi) = -Etotai- Now consider the result of a hypothetical measurement of the 
energy of one of the constituents. In equilibrium, the state of each constituent should be such 
that the average outcome of an energy measurement should be the same; that is, (Hi) = E, 
where E = N^ 1 E totai \. In other words, the equilibrium state of each constituent must lie on the 



energy surface £e = (V'l-^ilV'} = E} i n the pure-state manifold for that constituent. Since 

N is large, this will ensure that the uncertainty in the total energy of the composite system, as 
a fraction of the expectation of the total energy, is vanishingly small. 

It is convenient to describe the distribution of the various constituent pure states, on their 
respective energy surfaces, as if we were considering a probability measure on the energy 
surface £e of a single constituent. In reality, we have a large number of approximately 
independent constituents; but owing to the fact that the respective state spaces are isomorphic 
we can represent the behaviour of the aggregate system with the specification of a probability 
distribution on the energy surface of a single "representative" constituent. 

In thermal equilibrium the resulting distribution should be uniform over the energy surface 
£e since it must maximise the entropy. Therefore, the density of states is given by 

U(E) = J 5{H{ip) - E)dV T . (1) 

Here, V denotes the pure state manifold and dVr is the associated Fubini-Study volume element 
of r. Once Q(E) is specified, the entropy is given by S(E) = kB^n^(E). It follows that the 
temperature and the specific heat can be deduced from thermodynamic relations TdS = dE 
and C(T) = dE/dT. A short calculation shows that 

kBT= nWy and C( - T) = (2) 

The advantage of the present formulation over the traditional approach is that the entropy is a 
continuous function of the energy. As a consequence, thermodynamic functions such as those in 
([2]) are well defined for finite quantum systems. However, to justify the term "temperature" 
for the ratio Q/Q' we must show its properties are consistent with the requirements of 
thermodynamic equilibrium. For this purpose, consider two independent systems, each in 
equilibrium, with state densities [^(Ei)]^ 1 and [Q 2 (E 2 )] N2 • We let them interact for a period 
of time, during which energy e is exchanged. We then separate them and let them relax 
again to equilibrium. Because of the interaction the state densities of the systems are now 
[Qi(Ei + e/Ni)] Nl and [£l 2 (E 2 — £/N 2 )] N2 . The value of e is determined so that the total entropy 
S(E) = k B lnpl(Ei + e/N 1 )] Nl [n 2 (E 2 - e/N 2 )] N2 is maximised. This condition is satisfied if 
and only if e is such that the temperatures of the two systems defined according to ([2]) are equal. 
It follows that our definitions are thermodynamically consistent. 



3. Expressions for the density of states 

Let us now try to obtain a direct representation for the density of states £l(E) in terms of the 
energy eigenvalues. We consider first the two-level system with energy eigenvalues E\,E 2 . The 
Fubini-Study volume element for the pure state manifold is dVr = \ sin 9d9d<fi, where < 9 < tt 
and < <fi < 2tt. Since the energy expectation in a generic state = cos ^ 0|i^2)+sin 5 ^e'^l-E^) 
is E 2 cos 2 \ Q + Ei sin 2 \ Q = \{E 2 — E{)(1 + cos 0) + E\, the density of states is 

dA d^ r d9e- iX ( E{1+cosd ^- E ) S m9. (3) 
Jo Jo 

Here, we have made use of the integral representation for the delta-function, and we have also 
defined E = (E 2 — E\)/2. By use of the relation dAe _l5A A _1 sin(aA) = tt for a > b and = 
for a < b we thus deduce that Q(E) = tt/(E 2 — E{) for E\ < E < E 2 , and Q(E) = otherwise. 

An analogous calculation can be performed for a three level system. If we let Ei,E 2 ,E% 
denote the eigenvalues of the Hamiltonian, then the energy constraint for a generic state 



n(E) = SF L 



\ip) = sin g cos \ ip\E%) + sin | sin ^ yje^l-E^) + cos ^ ^ e 1T? | ) is given by E3 sin 2 ^ #cos 2 \ ip + 
E2 sin 2 ^ sin 2 ^ </? + E\ cos 2 ^0 = E. Since the Fubini-Study volume element in this case is 
dV r = ^ sin 0(1 — cos 0) sin</?d#d92cU;d?7, we carry out the relevant integration and obtain 

_ 7r 2 (£-£i) ^ ^ _ ir 2 (E-E 3 ) 

(E3 — E\)(E2 — .Ei) (£3 — Ei)(Es — E2) 

depending on E\ < E < E2 or E2 < E < E3. 

By pursuit of this line of argument we deduce more generally that the density of states £l(E) 
is given by a piecewise polynomial function of energy E. In particular, if the energy spectrum 
is nondegenerate, then we have the representation 

/ ^ n n+l n+l 

Here 1/a} denotes the indicator function = 1 if ^4 is true, and otherwise). To offer an 

intuition for the behaviour of the density of states, examples of £l(E) are shown in Figure [TJ 
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Figure 1. Density of states Q(E) given 
in (jHJ) associated with nondegenerate 
(n + l)-level systems for n = 3, 4, 5, 
with energy eigenvalues E^ = k (k = 
0, ... , n). For an (n + l)-level system, 
Q(E) is a piecewise polynomial of degree 
n — 1 , and is n — 2 times differentiable in 
energy E. 



Once the density of states £l(E) is obtained for the microcanonical equilibrium, thermal 
expectation values of physical observables in the corresponding canonical distribution can be 
computed by use of the canonical partition function Z(/3), which is the Laplace transform of 
£l(E). When energy eigenvalues are nondegenerate, we have 

n+l n+l 

k=i i=i^k' y 1 k ' 

A line of argument in Khinchin pQ for classical systems can then be applied here in the quantum 
context to prove that the parameter (3 appearing in ([6]) for the canonical partition function 
agrees with the microcanonical definition of temperature in (|2j). 



4. Quantum phase transitions 

An interesting consequence of the microcanonical framework, whether classical or quantum, 
is that the density of states in general need not be an analytic function for finite systems. 
In contrast, the partition function in the canonical counterpart is necessarily analytic. In 
other words, while it is necessary in the canonical framework to take thermodynamic limit 
to describe phase transitions, in the microcanonical formalism this is not the case. Therefore, 
an approach based on microcanonical equilibrium might provide an adequate description of the 
phase transitions for small systems observed in the laboratory. We note in this connection 
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Figure 2. A gas of quantum Ising chains. The gas consists of a large number of weakly 
interacting quantum molecules. Each molecule is modelled by a quantum Ising chain of three 
strongly interacting spin-i particles. 



that there are many classical systems for which finite-size phase transitions are predicted in 
microcanonical equilibrium • 

In quantum microcanonical equilibrium, the breakdown of analyticity of Q(E) gives rise to 
phase transitions in the sense that discontinuities in the higher-order derivatives of Q(E) emerge. 
Specifically, if we solve the first equation in ^) for the energy to obtain E(T), then for a system 
with n + 1 nondegenerate energy eigenvalues, the (n — l)-th derivative of the energy with respect 
to the temperature has a discontinuity. 

As an illustration we consider the specific heat C(T) for a gas of weakly interacting molecules, 
where each molecule is modelled by a strongly interacting chain of three Ising-type spins (see 
Figure [2]). The molecular Hamiltonian is H = —JY^k=i a z 'z +1 ~ ^ULi^: where o k z is the 
third Pauli matrix for spin k, and J, B are constants. For this system, the specific heat grows 
rapidly in the vicinity of the critical point T c = (2 J + B)/3ks, where the system exhibits a 
discontinuity in the second derivative of the specific heat. The plot of the specific heat is shown 
in Figure EJ along with the corresponding plot for a simple four-level molecular gas; the latter 
exhibits a second-order phase transition at the critical temperature k^T c = \e and critical 
energy E c = e, where e is the spacing of energy eigenvalues. 
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Figure 3. Specific heat for a nonde- 
generate four- level system (dotted line, 
n = 3, Ej = 0,1,2,3), and a quan- 
tum Ising chain having a four distinct 
degenerate eigenvalues (solid line, J = 
1/4, B = 1). The plots illustrate the 
existence of second-order phase transi- 
tions. In the quantum Ising chain, we 
have C(T) ~ (T - T c y 2 away from T c , 
whereas in the vicinity of T c we have 
C(T) ~ (T — T c )~ 13 for T > T c . 



5. Towards quantum grand microcanonical equilibrium 

In the foregoing discussion we have made use of the energy conservation property of the unitary 
evolution to introduce a quantum microcanonical hypothesis which asserts that in equilibrium, 



every quantum state with given energy E is realised with an equal probability. This hypothesis 
can be refined in the following manner, leading to what might appropriately be called the 
quantum grand microcanonical hypothesis. 

For a given quantum mechanical system there are n linearly independent conserved 
observables, where n + 1 is the Hilbert space dimensionality. Therefore, when an isolated 
quantum system with a generic Hamiltonian evolves unitarily, the associated dynamics exhibit 
ergodicity on the toroidal subspace T n C £e of the energy surface determined by simultaneously 
fixing the expectation values of the commuting family of observables (cf. [HI])- A theorem of 
Birkhoff [lj applies to show that the dynamical average of an observable can be replaced by the 
ensemble average with respect to a uniform distribution over T n . The density of states is then 
determined by the weighted volume of the subspace T n of the quantum state space. 

In the case of the energy observable, the conjugate variable is given by the inverse 
temperature. For other observables belonging to the commuting family, the associated conjugate 
variables can be thought of as generalised chemical potentials. In this respect, a refinement of 
the microcanonical postulate leads to an ensemble of the grand canonical form. 

Let us consider the simplest nontrivial example n = 2. We choose the two projection operators 
III = l-EaX-^il and II2 = l-^K-Ebl for the independent pair of commuting observables. Since 
these observables are conserved, we let the two constraints be (111) = p and (H2} = q. It follows 
from the resolution of identity that (II3} = 1—p — q. In terms of the usual parametrisation, these 
constraints read cos 2 \ = p and sin 2 ^ #sin 2 | ip = q, respectively. The generalised density of 
states Q(p, q) = f r ^((TIi) — p)5((Il2} — q)dVr can then be calculated to yield 

«(P, Q) = (T(p) - T(p + q - 1)) (T(q) - T(q - 1)) , (7) 

where T(x) = —1 for x < and T(x) = 1 for x > 0. We have Q(p, q) = tt 2 in the range 
< p, q < 1 and < 1 — p — q < 1. To establish its relation with the density of states $l(E) 
we solve the energy constraint pE\ + qE2 + (1 — p — q)Es = E for, say, p, then substitute the 
result in q), and integrate over q from to 1. The temperature of the system can then be 
obtained by differentiation. It would be of interest to further investigate properties of the grand 
microcanonical equilibrium for general systems, which in our view holds the promise for many 
new rigorous results in quantum statistical mechanics. 
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